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Abstract — This study consists of evaluating the performance 
and the accuracy of three filter bank design methodologies for 
the identification of the complex textural primitives of a SAR 
satellite image (4.6 meter resolution). The development of these 
filters is made in order to provide an accurate analysis of the 
available SAR images. For this purposes, we have used the 
following orthogonal spatial transformations: Karhunen Loeve 
Transform (KLT), Gabor Filters (GF) and Generalized 
Gaussian Modeling (GGM) by Wavelet Transform (WT), which 
allow the reduction of the correlation between data. The filter 
banks are represented by a combination of the optimal filters, 
selected for each texture according to the used transform. This 
paper shows the difficulty to choice the adequate design of the 
filter banks that allows an optimum identification of the textures 
of SAR images and classification of these images based on their 
filtered attributes. The results showed that the conventional 
filters used in the WT are not accurate for the heterogeneous 
textural identification, while high identification rate exceeding 
99% can be obtained with an adequate combination of the 
filters. The best filter bank is the one that allows the 
identification of all textures in spite of their dimension. 

Index Terms — Filter Banks, Image Textures, Texture 
Extraction, Filter Combination 

I. INTRODUCTION 

The identification of textures of satellite images has 
become a common practice used in many fields such as 
agriculture and forestry, and several works have been 
published on this topic. However, this task becomes harder 
when the textures are heterogeneous and irregular, requiring 
the use of appropriate methodologies and sometimes the 
adjustment of certain data. The Gabor filters, which are 
known for their robustness, have been widely used to a broad 
range of image processing tasks. 

The Gabor filters are very popular due to their conjoint 
space-frequency resolution which is optimal [8]. However, 
the design of the specific filters must be realized with high 
precision if the detection of the textural forms is required. 
Recently, Zand and al. [16] have used the Gabor and Curvelet 
approaches for the identification of the different textures of 
satellite images and modeling. However, as the Gabor and 
Curvelet transforms are designed for the regular shapes, only, 
a fitting method was applied to encode the information of the 
sub-bands in the polynomial coefficients, giving a texture 
feature vector with a maximum power of discrimination. For 
all these reasons and others, other approaches are proposed. 
Indeed, it was shown that it is possible to use the GGM of the 

Marrakchi Olfa Charfi, National Institute of Applied Science and 
Technologies, Member of GreenTeam-INAT LR17AGR01 Laboratory, 
Tunis, Tunisia 

Mbainaibeye Jerome, University of Doba, Doba, Chad 


local variances distributions of the energies of the sub-images 
issued from the WT to have a good classification of the 
different textures present within the remotely sensing images, 
amongst the SAR image [15]. With this method, it is also 
possible to develop the textures filter banks as we did in the 
following sections. 

The orthogonal transformations such as: KLT [3], WT [4], 
frames representation for textures [5], GF [9], were found to 
be more appropriate for classifying complex textures, 
allowing an optimal segmentation of the remote sensing 
images. It was observed that some descriptors such as the 
textural signatures of the SAR image textures, which are 
determined in the spatial representations using the KLT, WT 
and GF tools, are more efficient in characterizing such image 
textures [1]. Indeed, the KLT contributes to increase the 
capacity of the discrimination procedure of the objects inside 
the images. The information obtained from the new KLT axis 
constitutes the statistic information, which allows an easy 
identification of the studied texture. This image may be 
decomposed in the new orthogonal basis by filtering. The data 
obtained in the new basis are uncorrelated [2] [3]. 

The WT consists in analyzing a signal by using a localized 
centered function which may be translated and dilated [4], 
providing uncorrelated data [5] and allowing a 
space-frequency analysis identification of the global and the 
local textural informations. The WT is based on two 
conjugate filters: low pass filter and high pass filter, which 
are, also called quadrature mirror filters (QMF)[6]. With the 
WT tool, it is possible to have a decomposition of an image 
into sub bands called low pass sub band and high pass sub 
bands, which contains the details of the image [7]. Thus, the 
complexity of the data, their dimension and their inter 
relationships strongly influence the credibility of the results 
that can be obtained and define the degree of difficulty in 
analyzing the heterogeneous textures. In such analysis, it is 
important to localize the information at both, the spatial 
domain to characterize the orientations of the primitives and 
at the frequential domain to obtain the grey scales variation 
that can be used as descriptors; in this context, the best 
conjoint localization (spatial and frequential domains) may be 
obtained by the Gabor filtering [9] [10], with the cosine filter, 
which has a certain direction modulated by a Gaussian 
window [11][12]. 

In this paper, we are interested in the development of 
texture filter banks which are used to analyze the textures 
primitives and their attribute images. The remaining of this 
paper is organized as follows: In the section 2, we present the 
methods and materials used; the section 3 presents the results 
obtained; the section 4 presents the discussions of these 
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results; finally, the conclusion of this work is presented in the 

section 5. of the filters h tfc) and ,g(k) are calculated by the resolution 


II. Methods And Material 

A. Filter bank design based on transform 
A.l The KLT [1] 

The KLT image based transform is operated locally by 
considering a narrow neighboring, which is represented by a 
random vector X (its values are taken in the observation 
window where the dimension is k b x k b with k b x k b =N. The 
different realizations of X are obtained by moving the 
observation window through the image. For each value of X 
and for a given spatial position, is attributed its image Y by a 
linear transform represented by equation (1). The vector Y 
represents the image of the vector X in an orthonormal basis (f> 
given by equation (2) and contains N eigenvector of the X 
covariance matrix. 

Y = <p T X (1) 

Where: 






Y = 

U 

X = 

U 


The orthonormal basis is given by equation 2: 


of an equation system depending on and 

expressed by equation (5). 

=\ f 2Y k h(k)^‘(.2x - k) 

Wu) = jCfc3^Czjc - k) 

All couples of filters (h, g) designed in several scales 
represent the filter banks that we call Generalized Gaussian 
Model Filter Bank (GGM-FB). Figure 1 represents the 
scheme of the wavelet decomposition in one scale. 



Fig.l. One scale image wavelet decomposition. 


tp = 


0i,L 


0i.iv 

0iv.iv. 


( 2 ) 


According to [2], it is possible to establish a 
correspondence between the linear transform of a vector X by 
an eigenvector and the dimensional convolutional product 
of the observation of the window X(i, j) and a Finite Impulse 
Response filter (FIR) (/> n (i, j) represented by equation (3). This 
filter is constituted of the eigenvector coefficients reordered 
according to their positions in the observation window X(i, j). 

Yiif) = X[i,j) * ([',/) (3) 

Where (e,_, 0 is given by equation (4). 






(4) 


The image described by the variable Y (i, j) is called image 
attribute or filtered channel. A filter bank is composed of the 


FIR filters (/> n (i, j) with 1 < n < N. 


A. 2 The GGM [1][13] 

To local variance distributions of each detail subband 
energies of the image obtained by the WT at several scales 

may be associated to a scale function model . The first 
derivative of this function represents the wavelet function 

[4]. From the Gaussian model values of the function 


[7] and its first derivative the impulse responses 


A.3 The Gabor filters [13] 

The approach developed in this work, is based on a 
multichannel filtering of the image in a set of images 
containing the spatial characteristics at different scales. 
Additionally, the original image is decomposed in different 
planes corresponding to different frequential channels. Each 
of these images contains the textural characteristics appearing 
in a narrow spatial frequency band and orientation. This 
approach has the advantage to exploit the spatial interactions 
between pixels of neighboring at different scales. After 
filtering, we have selected the filters which are able to offer 
the capacity for efficiently detecting of the texture motives. 
The best simultaneous localization in spatial and frequential 
domain may be operated by the Gabor filtering (see equation 
6) applying a cosine filter with a certain direction modulated 
by a Gaussian window [11][12]. The radar image used in this 
work has only the real discrete values. The real form of the 
Gabor filter is defined by equation (7), where and >o y are 
the Gaussian standard deviations, f x and / v are the central 
frequencies, respectively in the x and y axis. To obtain 
different orientations, a rotation of the coordinate axis is 
made by using equation (8). 

& ^ 0x p{-£s-vti ex ? * + 


( 6 ) 


(7) 


Where 

(^/) is given by equation (8): 

( x i x/cosGj sinGA/xx 

W \-5m0j CO sGjAyJ { } 
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With 

f = 0 for e = 0 
! l =£ 0 for i += 0 

The GF has four parameters: the filter orientation, the 
central frequency f 0 to be extracted, the standard deviations 
a x and a y . For each texture, the GF are determined and 
permits its identification depending of the following 
variables: (cr x = rr y ) and f . Also, for each texture, a 

judicious choice of an optimal filter is realized among this set 
of filters according to the imposed criteria to its textural 
identification. The GF Filter Bank (GF-FB) is in this case 
composed of the combination of the optimal filters associated 
to each of the textures. 

B. Texture descriptor 
B.l Signature used for the KLT method [1] 

The signature vector C KLT attributed to each texture is 
formed by the energy eigenvalues of the texture; this vector is 

given by equation (9), where A £ is the i th energy eigenvalue 

extracted from the correlation matrix of a given texture and N 
represents the number of the energy eigenvalues. 

CjfZT = C^l .■ -Ij.- — Jfif ) (9) 

The number of the filter banks may be limited to a value 
lower than N according to the energy ratio ER (see equation 
10). ER is determined from the eigenvalues associated to the 
eigenvectors of the covariance matrix of the observation 

vector X, with > — > X N and p < N 
, . I* * 

ERtoti=^- (10) 

B.2 Signature used for GF method 

The signature vector, C GF , attributed to each texture is 
composed of two components determined on the filtered 
images by all the filters of the optimal filter bank. These 
components are expressed in terms of energy E and LI 
norm V given respectively by equations (12) and (13). 
This vector is given by equation (11), where E# defines the 
energy of the GF coefficients of texture 7) after filtering by 
the GF of frequency f], V X j is the LI norm obtained from 
the output of GF of frequency^. If a bank of four filters is 
used, the characteristic vector is written as indicated in 
equation (14), where E a (respectively Vu) represents the 
energy (respectively norm LI) of the GF coefficients for 
the frequency f of the texture i. 


CjTG 

= fci-Vj) 

(ID 




E <-r- 

PQ 

(12) 


E? =1 E? Juyjl 

(13) 

V v 

PQ 

L2 ‘ 

3 ■■ Efr ^Ll # ^L2 ' 

-DU) (14) 


B.3 Signature used for the GGM method 

The signatures in terms of energy (E) norm LI (V) are 
determined in the different detail subbands HL, LH and HH 
for the different scales of image wavelet decomposition. The 
signature vector V WT . GGM is determined by equation 15, where 
E t is the energy signature and V t is the LI norm signature; E x 
and Vi are determined in the i th detail subband at a given scale. 

Hi-T - S GM = ' ^LH ■' ' ■■ (15) 

C. Classification method 

The image classification based on minimum distance is 
realized according to the minimization of the mean square 
error (MSE) between the local statistics calculated in the 
scanning windows of k b x k b size on the filtered images and the 
textural signatures determined in the localized textural zones 
on the filtered images. The minimization of the MSE is given 
by equation (16), with Cj the j th component of the local statistic 

vector, Cj 1 the jth component of the characteristic vector of 

the texture 7) and n represents the number of the vector 
components which is a variable parameter depending on the 
selected transform method. 

min T .(EQM) = - cj l f (16) 

The image resulting from the minimum distance 
classification is an image of labels. The label represents the 
index i of the texture T x having a minimal MSE. 

The local statistics which are calculated in the scanning 
windows k b x k b are different from one transform to another: 

- For KLT, the local statistics are only represented by the 
average energy calculated according to equation (12) and 
determined in the observation windows on all the image 
attributes. In this case, the local statistic vectors c and the 
characteristics of the textures C have a number of components 
equal to the number of image attributes. The image is filtered 
by the filters of each filter bank (relative to a texture). All the 
image attributes by filter bank represent a list of images L x 
associated to a given texture T x . In the classification phase, the 
MSE is determined for each L x and for each position of the 
observation vector X. For a given position of X, the associated 
class label is those of the texture for which MSE X is the 
smallest for i = 1 to 4. This classification procedure is based 
on the method presented in [3]; 

- For WT-GGM, the local statistics are defined in terms of 
energy (see equation 12), and LI norm determined in the 
observation windows for all the images subbands for all the 
scales. In this case, the vectors (local statistics c and the 
textural characteristics C) have a number of components 
twice the number of the image detail subbands for K wavelet 
decomposition levels; 

- For GF, the local statistics are defined in terms of energy 
(see equation 12), and LI norm determined in the observation 
windows for all the filtered images by the filter banks. In this 
case, the vectors (local statistics c and the textural 
characteristics C) have a number of components twice the 
number of the filtered images. 

D. Evaluation criteria 

The filter banks are evaluated according the classification 
result by minimum distance of the SAR image textures in 
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terms of identification rate given by equation (17), with c the 
number of classes, I s the label of pixel s and p (I s =I/T S = T) is 
the conditional probability that the pixel s is classified in class 
I as it is in the texture T. 

!r(T) = £f =1 (p(4 = i/T s = T)/p(T s = T/I s = 0) X 
100 

(17) 

The validation of the filter banks is performed according to 
the chart shown in figure 2. 

The image is filtered by all the filters according to the 
applied transform method. From the filtered images and from 
the ground truth data location on these filtered images, the 
textural features are determined for each texture used for the 
airborne SAR image classification. 

The classification criterion used is the Mean Squared Error 
(MSE) between the local descriptor, calculated from the 
scanning windows ( k b x k b ) on filtered images (or attributes) 
and their textural signatures. 

The signatures are different from one method to another 
(KLT, GGM and GF). For some methods, a restriction on the 
number of filters in a filter bank is realized for the reduction of 
the computation amount in the filtering phase but keeping the 
maximum information of the image. 

The evaluation of the information is realized in this case by 
the computation of the energy ratio in percent (ER) expressed 
by equation (10). 

The selected filter banks are those which allow obtaining 
the best classification results in terms of identification rate 
according to equation (17) in the case of GGM and GF 
transforms. However, for KLT, the optimal filter bank is 
selected after the verification of the work presented in [3]. 



Fig. 2. Filter bank validation flow chart 


E. Datasets 

The textures used in this work are extracted from an image 
obtained by an airborne SAR in C band, in horizontal 
polarization and the resolution is 4.6 meters (Figure 4). This 
image covers a semi-arid region located in southern Tunisia. 
According to a ground truth campaign, soils in this region are 
essentially covered with scattered vegetation (V), olive trees 
(O.V), bare soils (BG.V) and urban areas (B.V). The thematic 
studies done on the zone have enabled us to identify the 
different textures on the airborne SAR image. This image has 
been the subject of several years of study since 1994. All 
texture classification operations have been performed on this 
image using several methods [13] [14]. Moreover, a 
comparative study of classification results with those of 
another radar image, acquired with RADARSAT satellite 
with lower resolution (12.5 meters), was performed before 
seeing the influence of the spatial resolution and its impact on 
the classification results [1]. We have deduced that the 
primitives in textures of the airborne SAR image are easily 
identified compared to those of the RADARSAT image. 
Hence, the interest is shown in these textures for developing 
and studying the performance of filter banks. 

The SAR image textures are irregular and consist of a 
common textural element which explains the strong 
correlation between the textures. For example, the textures 
O.V and BG.V have the “bare ground” as a common 
structural element; this is observed in figure 3. 



Fig. 3. Example of SAR image windows, (a) Olives on the 
ground covered with vegetation O.V (37 x 47pixels), (b) Built 
on ground covered with vegetation B.V(23 x 57pixels), (c) 
bare ground with the presence of vegetation BG.V ( 15 x 19 
pixels), (d) vegetation V ( 21 x 37pixels). 

The SAR textures are complex and are characterized by 
information redundancy. To reduce this textural correlation, 
the KLT [1], the WT and the GF [13] are used to represent 
these data in other spaces by conserving the global 
information of these textures [1] [14] [15]. 

III. Results 

A. Textural filter banks elaboration by the KLT 

For each texture, the filter banks and the associated textural 
signature vector (equation 9) are determined. The SAR image 
is filtered by the different elaborated filter banks and then, it is 
classified by the minimum distance method. Three parameters 
are adjusted for the optimization of the results: 

- The scanning window width k b x k b , used in the image 
classification phase; 

- The window of dimension l x l represented by a given 
textural texton used in the verification of the [3]; 

- The value attributed to energy ration ER. 
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The best SAR image classification results are obtained associated to the filters (or eigenvectors) are illustrated in 
for k b = 9,1 = 7 and ER = 90% figure 6 and they represent the textural signatures. 


.The best average identification rate is 80.2% as shown in 


Tab.2. 
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SAR innate 


Fig. 4. Airborne SAR image in C band, in horizontal 
polarization with resolution of 4.6 meters 


Fig.6. Variation of signature vector coefficients of SAR image 
texture for 1=7 and ER=90%. KLT method 



c d 

Fig.5. The first image attribute associated to the first 
more important energies for each texture (KLT method): 

(a) for O.V texture, (b) for B. V texture, (c) for BG. V 
texture and (d) for V texture 

Figure 5 illustrates the results of the KLT applied on SAR 
images and shows the first image attribute for each texture. 

This study shows that the number of filters varies from one 
texture to another (see Tab.l) and the computation time may 
be considered as a constraining parameter because the 
number of filtering operations is high. Only the eigenvalues 


Table. 1. Number of filters by texture for the SAR image 
transformed by KLT. ER = 90% and / = 7 


Textures 

O.V 

B.V BG.V V 

Number of filters 

15 

26 21 23 

Table.2 Identification ratio of SAR image textures in %. 

O.V B.V BG.V 

V 

Average 

Identification ratio (%) 

58.80 92.90 86.40 

82.70 

80.20 


B. Textural filter banks elaboration by the GGM 

The best classification results by minimum distance using 
GGM of the distributions of local variances in the detail 
subbands of SAR image was achieved for a size of scanning 
window k b = 11 and for two scales [15]. 

The scale function and wavelet function have been 
modeled on each scale of WT. By using these parameters and 
these models, and by solving the system of equations, the 
filters h and g were determined on each scale of WT for each 
of the textures in the image. Thus, for each texture are 
associated two filter banks: 1 bank by 2 scales of WT. Each 
filter bank is composed of three couples ( h,g ) extracted from 
the detail subbands: HL, LH and HH. Subsequently, the filter 
banks elaborated are used to realize the WT of the original 
image on two scales (Figures 7 and 8). 

The principle of the WT has been applied on the SAR 
image by using the filter banks h and g realized (see Table 3). 
Thus, to get the horizontal detail subband image (HL), we 
applied the high pass filter g HL , on the rows and then the low 
pass filter, h HL , on the columns of the original image (or the 
approximation subband if it concerns the second scale of 
WT). To get the vertical detail subband image (LH), we 
applied the low pass filter, h DV , on the rows and then the 
high-pass filter, g DV , on the columns of the original image (or 
the approximation subband if it concerns the second scale of 
WT) and finally, to get the diagonal detail subband image 
(HH), we applied the high pass filter, g DD , on the rows. This 
procedure is repeated according to the scale number in WT 
decomposition. 
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The best identification rates of the SAR image textures 
were obtained using the signatures of textures based on the 
detail subbands after application of WT using the filter banks 
of the texture B.V on one scale and for a scanning window 
size which is k b = ll(see Fig.9). Thus, the average rate of 
good classification of the four textures of SAR image reached 
91.31% (see Tab.4). Note that the identification rate obtained 
for BG.V texture is 100%. 



Fig.7. Example of scale functions obtained from the GGM 
of HH detail subbands of SAR image textures, (a) Texture 
O.V, subband HH, 1 scale, (b) Texture B.V, subband HH, 
1 scale, (c) Texture BG.V, subband HH, 1 scale, (d) 
Texture V, subband HH, 1 scale. 



Fig.8. Example of wavelet functions obtained from the 
GGM of HH detail subbands of SAR image textures, (a) 
Texture O.V, subband HH, 2 scales, (b) Texture B.V, 
subband HH, 2 scales, (c) Texture BG.V, subband HH, 2 
scales, (d) Texture V, subband HH, 2 scales. 


Table.3. Filter banks coefficients for the texture B.V of 
SAR image (one scale), (a) Filter h , (b) filter g 


Coefficients 
of filter h 

DH 

DV 

DD 

hd) 

-0.5234 

-0.5086 

-0.5883 

h(2) 

0.5593 

0.5633 

0.6716 

h(3) 

0.6207 

0.5766 

0.4871 

h(4) 

0.6482 

0.6648 

0.8396 

h(5) 

-0.5993 

-0.5918 

-0.7125 


Coefficients 
of filter g 

DH 

DV 

DD 

gd) 

0.0098 

0.0094 

0.0098 

g(2) 

-0.0091 

-0.0100 

-0.0124 

g(3) 

0.0144 

0.0145 

0.0150 

g(4) 

0.0003 

0.0017 

0.0068 

g(5) 

-0.0154 

-0.0156 

-0.0191 


Tab.4. Identification average rate of the SAR image 
texture for BF-MGG. 

k b = 11, B.V texture filter bank for one scale 


O.V 

B.V 

BG.V 

V 

Average rate (%) 

88.10 

82.50 

100 

94.70 

91.30 


C. Textural filter banks elaboration by EG 

The elaboration of Gabor filter banks requires two steps: 
firstly, we must determine, for each texture, different filters 
for its identification and furthermore for the enhancement of 
its textural patterns without alteration of their visual structure 
(see Tab.5). Then, several combinations of filters are realized: 
we choose a filter by texture, and then the selected filters are 
grouped to make a combination of filters. For each texture is 
attributed a textural signature (E, standard LI), determined on 
the image filtered through the optimal filter associated to this 
texture. The classification by minimum distance will indicate 
us what among the combinations of textural filters will 
achieve the best texture identification rate (see Tab.6). This 
method is more complex and requires a big computation time. 
It also requires a visual validation, (see Tab.6), so its 
realization cannot be automatic. 



Fig. 9. Variation of the texture identification rates of SAR 
textures by the differents filter banks: G-EB, GGM-EB 
and KLT-BF 


Table.5. The Gabor filter parameters chosen for each 
texture 


Textures 

ft values 

O.V 

fr={0.0 1 ,0.03,0.0 9 ,0. 11 ,0.2.0 3,0.4,0.5 } 
a = 0.5, 0=0 

B.V 

f<j={0.01 fl .09,02,0.3,0.4,0.5} 

I? = 0.5, 0=0 

BG.V 

fe= {0.01,0.09,0.11,0.2,0.3,0.5} 
o - 0.5, 6 - 0 

V 

fc={0.01.0.09.02,0.3.0.4.0.5} 
o = 0.5. 6 = 0 
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Table.6. Optimal textural identification rate of SAR image 
Gabor filter bank. 

frequencies = [0.09,0.09,0.01.0.09, 0 = 0 = , a = 0.5, 
k b =9. 


o.v 

B.V 

BG.V 

V 

Average 
identification 
rate (%) 

99.70 

98.60 

100 

99.60 

99.50 


Table 7. Variation of the scanning window size in 
function of the method used for an optimal 
classification result of SAR image 


Method 

KLT-FB 

GGM-FB 

G-FB 

h 

7 

11 

9 



Fig.10. The best texture classification results of SAR image, 
(a) KLT-FB, (b) GGM-FB, (c) G-FB. 

IV. Discussions 

The three methods used for the development of the filter 
banks are quite complex in their implementation. GF method 
applied to the SAR image is estimated to be faster than the 


GGM and KLT methods. In addition, it allows identifying all 
the textures of this image with classification rates up to 90% 
(see Fig.9).The complexity of the TKL method is due to the 
management of a very large number of filters by texture. 
However, the complexity of the GGM method is due the 
realization of the filters. 

The textural identification results show that the average 
rates of the SAR textures identification are greater than 85% 
(see Tab.2, Tab.4, Tab.6). Also, it has been shown that the 
nature of the texture affects the performance of the different 
texture identification methods (see Fig.10). Indeed, the size of 
the scanning window k b , representative of the texture, varies 
from an identification method to another (see Tab.7). 

V. Conclusion 

In this paper, different methods for the elaboration of the 
textural filter banks are presented. Only three methods among 
others are tested. The image used has been acquired with an 
airborne SAR sensor with high spatial resolution (4.6 m). 
Thus, the textures of this image are complex and present a 
mixture of primitives and therefore they are highly correlated. 
To reduce the degree of correlation, different orthogonal 
transforms were applied on our image by using the following 
methods: KLT, GF and the WT. We have shown that the built 
GF-FB permit us to reach an average textural identification 
rate of 99% of good classification, in the case of our SAR 
image textures. These results are more efficient than those 
obtained by the use of one of the transformations WT, FG or 
TKL methods without filtering operation. We have also 
mentioned that the effect of the size of the scanning window 
which is representative of the texture can gives the best 
classification results; however, the effect is variable from one 
method to another. This can be explained by the spatial 
transformation of the image. 
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